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ABSTRACT 

We have tested the ability of driven turbulence to generate magnetic field structure from a 
weak uniform field using three dimensional numerical simulations of incompressible turbulence. 
We used a pseudo-spectral code with a numerical resolution of up to 144 3 collocation points. We 
find that the magnetic fields are amplified through field line stretching at a rate proportional to the 
difference between the velocity and the magnetic field strength times a constant. Equipartition 
between the kinetic and magnetic energy densities occurs at a scale somewhat smaller than the 
kinetic energy peak. Above the equipartition scale the velocity structure is, as expected, nearly 
isotropic. The magnetic field structure at these scales is uncertain, but the field correlation 
function is very weak. At the equipartition scale the magnetic fields show only a moderate 
degree of anisotropy, so that the typical radius of curvature of field lines is comparable to the 
typical perpendicular scale for field reversal. In other words, there are few field reversals within 
eddies at the equipartition scale, and no fine-grained series of reversals at smaller scales. At 
scales below the equipartition scale, both velocity and magnetic structures are anisotropic; the 
eddies are stretched along the local magnetic field lines, and the magnetic energy dominates the 
kinetic energy on the same scale by a factor which increases at higher wavenumbers. We do not 
show a scale-free inertial range, but the power spectra are a function of resolution and/or the 
imposed viscosity and resistivity. Our results are consistent with the emergence of a scale-free 
inertial range at higher Reynolds numbers. 

Subject headings: ISM:general-MHD-turbulence 



1. INTRODUCTION 

There are at least two distinct types of mag- 
nctohydrodynamic (MHD) turbulence. When the 
external large scale magnetic field is strong, the 
resulting turbulence can be described as the non- 
linear interactions of Alfven waves (e.g., Goldre- 
ich and Sridhar 1995, 1997). In contrast, when 
the external field is weak, MHD turbulence near 
the scale of the largest energy containing eddies 
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will be more or less like ordinary hydrodynamic 
turbulence with a small magnetic back reaction. 
In this regime, the turbulent eddy turnover time 
at the large scale (L/V) is less than the Alfvenic 
time of the scale (L/B), where V and B are r.m.s. 
velocity and magnetic field strength (divided by 
(47T/9) 1 / 2 ) respectively and L is the scale of energy 
injection or the largest energy containing eddies. 

Various aspects of weak/zero external field 
MHD turbulence have been studied both theo- 
retically and numerically. One of the most im- 
portant issues in this regime is the generation of 
large scale fields. Since large scale magnetic fields 
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are observed in almost all astrophysical objects, 
the generation and maintenance of such fields is 
of great importance. In mean field dynamo the- 
ory (see Moffatt 1978; Parker 1979), turbulent 
motions at small scales are biased to create a 
non-zero electromotive force along the direction 
of the large-scale magnetic field. This effect (the 
'a-effect') works to amplify and maintain large- 
scale magnetic fields. Whether or not this effect 
actually works depends on the structure of the 
MHD turbulence, especially on the mobility of 
the field lines. For example, Vainshtein and Cat- 
taneo (1992) have argued that when equipartition 
between magnetic and kinetic energy densities oc- 
curs at any scale larger than the dissipation scale, 
the mobility of the field lines and the a effect 
will be greatly reduced. An alternative approach, 
which may avoid such difficulties, is to appeal to 
small scale field line stretching alone to generate 
strong magnetic fields (Batchelor 1950, Kazantsev 
1967). In this model large scale magnetic fields 
are actually composed of parallel fibrils of random 
polarity, i.e. the large scale field is defined by 
the preferred axis of the magnetic field direction, 
but there is little large scale flux in such fields. 
Whether or not the mean field dynamo mechanism 
works, it is interesting to contemplate whether or 
not such fine-grained fields can be produced in the 
absence of any a effect. 

In addition, there are at least two other points 
of contention regarding the nature of magnetic 
fields generated in MHD turbulence. First, the 
mobility of magnetic fields will be affected by their 
intermittency. Various models for the generation 
and dynamics of flux tubes have been proposed, 
motivated in part by observations of the solar pho- 
tosphere (e.g., Vishniac 1995a, 1995b; Branden- 
burg at al. 1995). However, so far numerical sim- 
ulations have tended to show only modest levels 
of intermittency. Second, there is the nature of 
the energy cascade and power spectrum of MHD 
turbulence. There are numerous analytical mod- 
els for dealing with MHD turbulence: the eddy- 
damped quasinormal Markovian (EDQNM) ap- 
proximation (Pouquet et al. 1976), the renormal- 
ization group technique (Fournicr at al. 1982), and 
others, including recent work by Goldreich and 
Sridhar (1995, 1997) in which they treated nonlin- 
ear eddy interactions as a series of interactions be- 
tween marginally nonlinear Alfven waves. The last 



is perhaps the most promising. The energy spectra 
of the MHD turbulence is one of the most poorly 
understood fields in astrophysics. Early work by 
Iroshnikov (1963) and Kraichnan (1965) obtained 
a k~ z / 2 spectrum for both magnetic energy and 
kinetic energy in the presence of a dynamically 
significant magnetic field. However, this work was 
based on assumption of isotropy in wavenumber 
space, which is difficult to justify unless the mag- 
netic field is very weak. Goldreich and Sridhar's 
model assumes a critical level of anisotropy, such 
that magnetic and hydrodynamic forces are com- 
parable, and predicts a A: -5 / 3 spectrum for strong 
external field turbulence (which should apply to 
very small scales within any MHD turbulent cas- 
cade). Solar wind observations, which are well 
within the strongly magnetized regime, show an 
energy spectrum E(k) <~ fc -170 (Lcamon et al. 
1998). Numerical studies have only recently been 
able to address this question. Recent work by 
Cho and Vishniac (2000) and Maron and Goldre- 
ich (1999) seems to support the scaling laws of 
Goldreich and Sridhar. 

Here we will concentrate on the structure of 
magnetic fields in MHD turbulence when the ex- 
ternal field is weak. In the case of hydrodynamical 
turbulence, the energy cascades to smaller scales. 
If we introduce a weak uniform magnetic field, 
turbulent motions will stretch the magnetic field 
lines and divert energy to the small scale magnetic 
field. As the magnetic field lines are stretched, the 
magnetic energy density increases rapidly, until 
the generation of small-scale magnetic structures 
is balanced by the magnetic back reaction at some 
scale between L and the dissipation scale. This 
will happen when the magnetic and kinetic energy 
densities associated with a scale I are compara- 
ble so that Lorentz forces resist further stretching 
at or below that scale 2 . However, stretching at 
scales larger than I is still possible, and the mag- 
netic energy density will continue to grow if I can 
increase 3 . Eventually, a final stationary state will 
be reached. 

This discussion raises two questions. What is 
the scale of energy energy equipartition? What 

2 The first energy equipartition scale I can be larger than the 
dissipation scale when the equipartition occurs before the 
energy cascade reaches the dissipation scale. 

3 Our results show that I does increase and the magnetic 
energy continues to grow. 
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is the magnetic field structure on this scale? The 
answer to the latter question will depend on the 
nature of diffusive processes acting on the mag- 
netic field. First, suppose that magnetic field lines 
are unable to smooth the tangled fields at small 
scales. Then, as a result of the turbulent energy 
cascade and the subsequent stretching of magnetic 
field lines, magnetic fields have thin fibril struc- 
tures with many polarity reversals within the en- 
ergy equipartition scale I. Consequently, magnetic 
structures on the equipartition scale are highly 
elongated along the local magnetic field direction. 
This is the kind of picture one obtains by con- 
sidering passive advection of magnetic fields in a 
chaotic flow (for a review see Ott 1998). The de- 
gree to which this can be applied to a realistic, 
highly conducting fluid is controversial. On the 
other hand, if we assume MHD turbulence is al- 
ways capable of relaxing tangled field lines at small 
scales, then we expect eddies at the final equipar- 
tition scale to be nearly isotropic. That is, the 
typical radius of curvature of field lines will be 
comparable to the scale of field reversal. This pic- 
ture docs not imply a clear expectation for the 
scale of equipartition in the stationary state. As 
we will see later, our simulations support the lat- 
ter picture but with the energy equipartition scale 
near, but somewhat smaller than, the largest en- 
ergy containing eddy scale. 

In this paper, we will interpret the results of 
numerical simulations of MHD turbulence with a 
weak external field in terms of these models. In §2, 
we discuss the numerical method we used and per- 
formance of the code. In §3, we present the results 
of simulations. Finally, §4 contains a discussion of 
the results and our conclusions. 

2. NUMERICAL METHOD AND PER- 
FORMANCE OF THE CODE 

We used a pseudo-spectral code to solve the 
MHD equations in a periodic box of size 2n: 

— = (V x V) x V - (V x B) x B+^V 2 V+f +VP', 

(1) 

— = B VV- V- VB + ryV 2 B, (2) 

V-V = V-B = 0, (3) 

where f is a random forcing term with unit corre- 
lation time, P' = P + V • V/2, V is the velocity, 



and B is magnetic field divided by (A-Kp) 1 / 2 . Thus 
the field B is, in fact, the Alfvenic velocity. The 
Alfvcn velocity of the background field, B , is set 
to 10 -3 ' 5 in all simulations. Throughout the pa- 
per, we consider only cases where viscosity is equal 
to magnetic diffusivity: 

v = r). (4) 

In the pseudo-spectral approach, the actual calcu- 
lations are performed in wavevector space. The 
nonlinear terms are evaluated in real space using 
Fourier transformed variables and their derivatives 
and then transformed back into their Fourier com- 
ponents. The forcing term consists of 21 Fourier 
components with 2 < k < \[V1. The peak of 
energy injection is at k s=a 2.5(= k^). We ad- 
justed the amplitudes of the forcing components 
so that V ~ 1. We use exactly the same forcing 
terms for all simulations. Therefore, the external 
magnetic fields are very weak for all simulations 
(i.e. B <C V). Each forcing component con- 
sists of two parts: a linearly polarized component 
and a small circularly polarized component with 
a preferred helicity. The latter was introduced to 
provide non-zero helicity injection to the turbu- 
lence. The resulting fluid helicity is 0.3 ~ 0.4 for 
all runs. We use an appropriate projection oper- 
ator to calculate the VP' term in Fourier space 
and also to enforce the divergence-free condition 
(V • V = V • B = 0). We use up to 144 3 collo- 
cation points. We use the integration factor tech- 
nique for kinetic and magnetic dissipation terms 
and the leap-frog method for the nonlinear terms. 
At t = 0, the magnetic field is uniform and the 
velocity is spread between 2 < k < 4 in wavevec- 
tor space. We give the parameters for each run in 
Table 1. 

Either physical viscosity (and diffusivity) or hy- 
perviscosity (and hypcrdiffusivity) is used in the 
dissipation terms (see Table 1). The power of hy- 
perviscosity is set to 8, such that the dissipation 
term in the above equation is replaced with 

-^(V 2 ) 8 V, (5) 

where f§ is determined from the condition iSh(N/2) 2h At « 
0.5 (see Borue and Orszag 1995,1996). Here At 
is the time step and N is the number of grids in 
each direction. The same expression is used for 
the magnetic dissipation term. 
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When f = v = 7] = 0, that is, without forcing or 
dissipation, the total energy drops by 0.87 percent 
after 13 time units, apparently due to truncation 
errors. In the dissipative case the injection and 
dissipation of energy are well balanced after the 
magnetic energy reaches the stationary state: 

< e >=< D K > + < D M >, (6) 

where e is the energy injection rate and Dk and 
Dm are the kinetic and the magnetic energy dis- 
sipation rates respectively. The angle bracket 
stands for an appropriate space-time average, 
which is taken after the magnetic energy reaches 
the stationary state. 

Before we describe the third test, wc define sev- 
eral transfer functions. Imagine velocity compo- 
nents in a thin (i.e., thickncss=l) spherical shell 
of radius k in Fourier space. The work done to 
the velocity components by magnetic fields can be 
written as follows: 

- ]T V(p)*.[(V^TxB](p), (7) 

fc<| P |<fc+i 

where hatted variables are Fourier space variables 
and '*' denotes a complex conjugate. This relation 
is derived from the second term on the right-hand 
side of equation (1). This quantity is the negative 
of the energy transferred to the magnetic fields 
from the velocity components per unit time: 

W*0 = + E V(p)* • [(V >TbTx B](p).. 
fe<lpl<fc+l 

(8) 

Note that the subscripts of T x \y read Ho X from 
K' Similarly, from the first term in the right-hand 
side of equation (1), we obtain T v \ Vk (k): 

T v \ Vk (k) = - £ V(p)*-[(VrvTxV](p). 

fe<|p|<fe+i 

(9) 

Now, consider the magnetic components in the 
same shell. From the first term in the right-hand 
side of equation (2), we obtain the energy trans- 
ferred to the velocity fields from the magnetic com- 
ponents per unit time: 

T v \ bk (k) = - £ B(p)*.[B^VV](p). (10) 

fe<| P |<fe+i 

the work done to the magnetic components by ve- 
locity fields: Again, this quantity is the negative 



of the work done to the magnetic components by 
velocity fields per unit time. 

In our third test, we checked the energy budget 
of the magnetic fields. The net energy transferred 
to magnetic fields from velocity fields can be calcu- 
lated from either T b \ Vk (k) or T v \ bk (k). This energy 
will disappear through diffusive damping. There- 
fore, the numerical scheme should satisfy 

Tb\ v = —T v \b = D M , (11) 

where 

n lv =Y^T blVk (k), (12) 

fe=0 

T v \ b = E T v\b k ( k )> ( 13 ) 
fe=0 

and k max is the largest wave number in Fourier 
space and is equal to N/2. Fig. 1 shows that our 
code satisfies equation (11). The two curves Tb\ v (t) 
and — T v ib(t) coincide exactly. It is interesting to 
note that there is a time delay between T b \ v (= 
T v \b) and Dm, which is a measure of the energy 
cascade time. 

3. RESULTS 

3.1. Generation of Magnetic Fields 

We list the results of our simulations in Table 
1. We obtained V 2 , B 2 , e, and Dm by averag- 
ing over (ii,i2) using all available data. However, 
T v]v2 (= T v]vk (k = 2)) and T b[v2 (= T b \ Vk {k = 2)) 
were calculated from a sparse subset of the data. It 
is important to note that, unless otherwise noted, 
these time averages arc taken after the turbu- 
lence has reached the stationary state. Conclu- 
sions based on these averages do not apply to the 
initial growth phase of the magnetic field. 

Fig. 2 shows time evolution of kinetic and mag- 
netic energy density. All the simulations have sim- 
ilar kinetic energy densities. However, the mag- 
netic energy density obviously depends on the 
ohmic diffusivity n. After an initial growth phase, 
the magnetic energy reaches a stationary state. 
The initial growth rate of the magnetic energy de- 
pends on diffusivity. Even though it is not clearly 
shown in the figure, the initial growth phase con- 
sists of two stages for Run 144A, the hypervis- 
cosity run. In the first stage, which begins at 
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t=0 and ends after a few dynamical times, the 
growth rate is very fast. At this stage the mag- 
netic energy grows through the stretching of mag- 
netic field lines. Growth is fast because there is 
no significant back reaction by the magnetic fields. 
Fig. 3 shows that stretching is initially most ac- 
tive near, or somewhat larger than, the dissipa- 
tion scale (cut-off scale) and the magnetic energy 
spectrum peaks at this scale. As the magnetic en- 
ergy grows, the magnetic back reaction becomes 
important at the dissipation scale. When energy 
equipartition is reached at this scale, the stretch- 
ing rate slows down and a second stage of slower 
growth begins. Fig. 3 shows that during this stage 
the peak of the magnetic power spectrum moves 
to larger scales. Fig. 2 shows the second stage 
ends at t = 30 ~ 40 for Run 144A. Other runs 
with physical viscosities (Run 128A, 96A, 64A, 
etc.) also show similar behavior. However, the 
detailed evolution seems different from the hyper- 
viscosity case. For example, the diffusivity also 
affects the net growth rate in the second stage. 

In Fig. 4a, we plot the kinetic and magnetic 
energy densities as functions of v (= rj). Error 
bars are estimated from 

2<* 2 >-<* >2 , (14) 

V (*2 — t\)/t corr 

where < • > denotes the time average and X = V 
and B. As noted above, this average is taken af- 
ter the turbulence has reached its stationary state. 
The correlation time t corr ~ 3 is calculated di- 
rectly. We see that the kinetic energy is almost 
independent of v (= r/). The constancy of the 
kinetic energy density is the result of the mag- 
netic field. In purely hydrodynamic simulations 
the kinetic energy density increases as v decreases 
(see REF2 and REF3 in Table 1). On the other 
hand, the magnetic energy increases as v (= 77) 
decreases. This figure suggests that at small r\ 
even a weak diffuse magnetic field can lead to a 
strong final state. In the case of hyperviscosity, 
the magnetic energy is almost half of the kinetic 
energy at late times. This is the result of field 
line stretching, and not the alpha effect (see next 
paragraph and §4). The energy injection rates e 
(= f • V = Dk + Dm) are independent of v (= rj) 
and in the range 0.161 < e < 0.166 (Table 1). 

In Fig. 4b, we plot an empirical relation for 
the generation of magnetic fields in the stationary 



state: 

Dm ~ (t~ c t) b2 ' (15) 

where c is a constant (~ 0.63). It is a somewhat 
dangerous to assume that Run 144A, the hyper- 
viscosity run, is equivalent to a the case of vanish- 
ing viscosity. In particular, Fig. 5 shows a sharp 
change in the energy transfer to the magnetic field 
near the dissipation scale. This makes it difficult 
to extend this relationship to the case of vanishing 
diffusivity. 4 . 

What is the origin of this correlation? In the 
stationary state, the magnetic dissipation (Dm) 
is balanced by the net energy transferred to the 
magnetic fields from the velocity fields. The right- 
hand side therefore tells us that the net energy 
transferred to the magnetic fields is proportional 
to large-scale eddy turnover rate (V/L) minus 
an Alfvenic frequency (B/L) times a constant. 
The large-scale eddy turnover rate is equal to the 
stretching rate of the magnetic field when the back 
reaction is zero. We identify the second term on 
the right-hand side of this equation as the effect 
of the magnetic back reaction. 

This interpretation is supported by the spec- 
trum of the energy transfer rate, T b \ Vk (k). The 
function T b \ Vk (k) is the energy lost by velocity 
components within a unit shell of radius fe in 
wavevector space through interactions with the 
magnetic field. When this quantity is positive, 
it means that energy is being transferred from ve- 
locity components in the unit shell to the mag- 

4 In Run 144A, the magnetic dissipation is given by 
^2 fc 1 6 1 B 1 2 . Although we can calculate this quantity di- 
rectly, we should not use this for equation (15) because 
Dm in equation (15) actually represents the amount of en- 
ergy transferred to the magnetic field on large scales. If we 
use a physical viscosity (and diffusivity) , the two quantities 
are same. However, in hyperviscosity runs, the two quanti- 
ties diverge because nonlinear processes occurring near the 
dissipation cut-off are presumably unphysical. The dotted 
curve in Fig. 5 represents the amount of energy transferred 
to the magnetic field in Run 144A. The curve shows a non- 
negligible change after the dissipation cut-off (fc^ ~ 50), 
which means that nonlinear processes do occur within the 
dissipation range (i.e. fc > fc^). If we take the value of Iljj,, 
at k ~ fcd for Dm and substitute it into equation (15), then 
we obtain (y — B/l.6)B / Dm = 1-24, in good agreement 
with the other data shown in Fig. 4b. This suggests that 
equation (15) is is true even in the limit of vanishing viscos- 
ity (and diffusivity), but given the existence of anomalous 
energy transfer near the dissipation scale, this conclusion 
may be premature. 
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netic field. In Fig. 5a, we plot T v \ Vk (k), T b \ Vk (k) 
and T v \ bk {k) for Run 144A. From the spectrum of 
T b \ Vh {k) we can see that magnetic fields are driven 
by energy extracted from the large scale velocity 
components (cf. Kida et al., 1993). We can ex- 
plain this result if wc suppose that the magnetic 
field energy is generated by field line stretching 
at large scales. On the other hand, the magnetic 
fields have a small rate of net energy transfer to 
small scale motions, implying that on small scales 
the turbulent motions can be described as nonlin- 
ear Alfven waves, with a rough balance between 
kinetic and magnetic energy densities. 

In Fig. 5b, we plot the spectra of energy fluxes. 
The energy flux H b i Vk (k) is the amount of en- 
ergy transferred from velocity components whose 
wavenumbers are less than or equal to k to the 
magnetic components of all wavenumbers: 

IW*0= / T blVp (p)dp. (16) 
Jo 

We define other fluxes similarly. The graph of 
H b \ Vk (k) rises rapidly at large scales (k < 10) 
and reaches the value of energy injection rate e 
(<~ 0.16) at k <~ 10. This means that the energy 
injected to drive the turbulence (~ 0.16 per unit 
time) is absorbed at large scales by the magnetic 
fields. On smaller scales (k > 10) the velocity 
field does not drive the magnetic fields, but instead 
absorbs energy. In summary, the rise of n b |„ fc (k) 
at small wavenumbers is due to field line stretch- 
ing and the decrease at large wave numbers is a 
sign of Alfvcnic turbulence at small scales. Conse- 
quently, at small scales, we can treat the turbulent 
motions as fully developed MHD turbulence, even 
when external fields are weak. 

Before proceeding, it is useful to take note of a 
property of hyperviscosity simulations. The quan- 
tity TL v \ Vk (k) + Hb\ Vk (k) is the total energy trans- 
ferred from large scale velocity fields (|k| < k) to 
either magnetic fields or small scale velocity fields. 
When there is no dissipation, the energy injected 
by the driving force is either transferred to mag- 
netic fields or transferred to small scale velocity 
fields: 

e = n, K (fc) + n bK (fc). (17) 

That is, the quantity Yl v \ Vk (k) + U b ^ Vk (k) has to 
be a constant for a hyperviscosity simulation. In 
fact, in the hyperviscosity simulation (Run 144A), 
this quantity is constant for k < 50. 



In Fig. 6, we plot T V \ V2 and T b \ V2 . Here T V \ V2 
(= T v \ Vk (k = 2)) is proportional to the rate at 
which the energy of large scale eddies is transferred 
to small scale velocity fields. On the other hand, 
T b \ V2 (= T b \ Vk (k — 2)) is proportional to the rate 
at which energy of large scale eddies is transferred 
to magnetic fields. Note that k = 2 corresponds to 
the scale of largest energy containing eddies. Fig. 
6a shows 

T V]V2 ~ V\ (18) 

which suggests that large scale kinetic energy (V 2 ) 
is transferred to smaller scales within an eddy 
turnover time (L/V). That is, turbulence at large 
scale is broadly similar to ordinary hydrodynamic 
turbulence and the classical energy cascade model 
seems to work. In the case of hyperviscosity, the 
figure suggests that the energy transfer rate may 
not exactly follow the scaling relation given above. 
However, as we can see in the figure, the deviation, 
if any, will be very small. 
Fig. 6b shows 

T b{v2 ~ B 2 V. (19) 

When magnetic field is generated through stretch- 
ing of field lines, T b i V2 will be proportional to the 
stretching rate at large scales. Although equation 
(15) implies that overall stretching rate is propor- 
tional to (V — const ■ B), Fig. 6b suggests that the 
stretching rate scales as B 2 V at the largest energy 
containing eddy scale. That is, it is proportional 
to the large eddy turnover rate V/L. 

3.2. The Structure of Turbulence 

We plot energy spectra in Fig. 7. For all simu- 
lations using physical viscosities (runs 128A, 96A, 
72A, and 64A-D) the kinetic energy spectra are 
almost independent of v at large scales (k < 10). 
However, this is not true in the case of hypervis- 
cosity (Run 144A). Magnetic energy spectra peak 
near, but somewhat smaller than, the energy in- 
jection scale. The position of the peaks depends 
on magnetic diffusivity. We see that the location 
of the peaks moves towards the largest scale as 
diffusivity increases. In the case of hyperviscos- 
ity (and hyperdiffusivity) simulation, both kinetic 
and magnetic spectra show the well known bot- 
tleneck effect (see, for example, Borue and Orszag 
1996) for k >^ 20. This effect is characterized 
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by a local enhancement of the energy spectra at 
scales somewhat larger than the dissipation cutoff. 

We define the energy equipartition wavenumber 
k e such that 

/ d P E M (p)= / dpE K (p). (20) 

Jo J fc c 

In this definition, energy equipartition between 
large scale magnetic fields and small scale velocity 
fields occurs at k = k e . (Therefore, k e is different 
from the scale of equipartition between small scale 
magnetic energy and small scale kinetic energy. A 
good measure for the equipartition scale between 
small scale energies is the wavenumber at which 
Ex{k) — EM{k). In the case of hyperviscosity, 
the two wavelengths are similar.) From Table 1, 
we can see that k e is not very sensitive to the value 
of viscosity. In all cases, the equipartition scale 
(<~ l/k e ) is somewhat smaller than the peak of 
the kinetic energy spectrum (1/fc^ ~ 1/2.5). The 
magnetic fields associated wavenumbers greater 
than k e act like a strong uniform external mag- 
netic field for eddies smaller than ~ l/fc e , and the 
structure of turbulence at smaller scales will be 
similar to Alfvenic turbulence. On the contrary, 
turbulence at scales larger than ~ l/fc e will be 
more or less like ordinary hydrodynamic turbu- 
lence. This argument is supported by the fact that 
the value of k e is very close to the peak of Tlb\v k (k) 
for all runs. For example, Fig. 5b shows n fc | t , fc (k) 
peaks at k <~ 10 for Run 144A. Table 1 shows 
that k e <~ 8.4 for the same run. (Our analysis also 
shows that, for Run 96A, k e <~ 6.7 and U b \ Vk (k) 
peaks at k ~ 8. And, for Run 64A, k e <~ 6.6 and 
Hf J i Vk (k) peaks at k ~ 8...) This suggests that, be- 
low the equipartition scale, the turbulent fluid is 
in a state of nonlinear Alfvenic turbulence, with a 
fixed ratio between the kinetic and magnetic ener- 
gies associated with perturbations on these scales. 

We see that a small external field leads to lo- 
cal magnetic energy density which, in the limit of 
very small diffusivity, is comparable to the kinetic 
energy density, with a typical scale only slightly 
smaller than the scale at which the turbulence is 
driven. This still leaves the question of whether 
the magnetic field is characterized by numerous 
polarity reversals within each eddy. We have ex- 
amined this question by considering the magnetic 
field second order structure function. A detailed 
description of our method can be found in Cho and 



Vishniac (1999). Here we note only that we define 
the structure function in cylindrical coordinates as 

F(p,s) = <<|Y(ri)-Y(r 2 )| 2 >, (21) 

where Y can be either the magnetic field or the 
velocity. The coordinates are defined relative to 
the local field direction, with p being the distance 
perpendicular to (B(ri) + B(r 2 ))/2 and z the dis- 
tance parallel to it. 

We plot the results for run 144A in Fig. 8. The 
shape of the contours represents the average shape 
of the eddies, which in turn reflects the degree of 
anisotropy along the local field direction. As ex- 
pected, the velocity field shows isotropy at large 
scales and anisotropy at small scales. The mag- 
netic fields show an insignificant amount of cor- 
relation at large scales and strong anisotropy at 
small scales. At intermediate scales, comparable 
to the equipartition scale, the magnetic field shows 
a factor of two difference in correlation lengths 
along local field line directions and across them. 
The implication is that the field is not composed 
of numerous parallel fibrils with frequent polarity 
reversals on small scales, as envisioned by Batche- 
lor (1950) and Kazantsev (1967). Instead, the field 
is largely smoothed by small scale turbulent diffu- 
sion, and shows only one or two polarity reversals 
at its typical scale of organization. 

4. DISCUSSION AND CONCLUSIONS 

We have presented results of MHD numerical 
simulations for unit magnetic Prandtl number (i.e. 
v = r\). A weak diffuse magnetic field was used as 
an initial condition, and we tested the sensitiv- 
ity of our results to dissipative effects by varying 
the viscosity and numerical resolution in our sim- 
ulations. We have found that magnetic fields are 
amplified through field line stretching at a rate 
proportional to (V — cB)/L, where c ~ 0.63. It is 
not certain whether or not the hyperviscosity run 
follows this relation. We have also shown that, 
in the limit of v (= rj) — > 0, the magnetic field 
reaches a final stationary state where the mag- 
netic energy density is comparable to the kinetic 
energy. Since the simulations using physically re- 
alistic functional forms for the diffusivity are very 
far away from this limit, we are not able to ex- 
trapolate to a definite final ratio of Eb/Ek, but 
it should be at least as large as the value reached in 
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run 144A, <~ 0.6. The typical scale of the magnetic 
field is slightly smaller than the scale of the largest 
energy containing eddies. At or above the equipar- 
tition scale, eddies are almost isotropic. However, 
eddies smaller than the equipartition scale show 
elongation along the local magnetic field lines. 

When we vary viscosity (and diffusivity) explic- 
itly, physical quantities scale as follows: 

1. V 2 ~ const 

2. D M ~ (V - cB)B 2 , c^constant 

3. T v \ Vk (k = 2) ~ V 3 

4. T blVk (k = 2)~VB 2 . 

The turbulent dynamo effect did not play a 
role in these simulations. There is no sign of 
the spontaneous generation of a large scale field, 
above what we would expect by extrapolating the 
magnetic energy power spectrum to wavenumbers 
smaller than the energy peak. Our results are 
most plausibly explained in terms of random field 
line stretching within eddies. We have tested this 
conclusion by conducting additional simulations 
with no imposed helicity. The results were almost 
identical to the simulations presented here. How- 
ever, this does not allow us to draw any general 
conclusions about the possibility of mean-field dy- 
namos in astrophysical objects. We note that the 
average helicity for all the runs listed in Table 1 
is in the range 0.3 to 0.4. Since the smallest al- 
lowable wavenumber is 1, this implies a maximum 
mean- field growth rate of order 0.1 times the eddy 
correlation time. By comparison, the turbulent 
dissipation rate is (V 2 )/3 (or <~ 0.3) times the 
eddy correlation time. The obvious conclusion is 
that no dynamo effect is expected. The energy 
injection scale here is too close to the size of the 
computational box. 

The isotropic structure of eddies at large scales 
in our simulations implies something about diffu- 
sive processes in MHD turbulence. In particular, 
it suggests the presence of an effective magnetic 
diffusivity which is tied to the large scale eddy 
size rather than the resistive scale. This is an 
old idea (see, for example, Parker 1955), and is 
equivalent to the notion that the usual turbulent 
diffusion coefficients can be substituted for ohmic 
resistivity. In its original form it was based on 
a picture of small scale field line mixing, which 



has been convincingly criticized by Parker (1992). 
However, it can also be justified by appealing to 
rapid reconnection, that is, field line reconnection 
at rates which are comparable to an eddy turn over 
rate and much more rapid than estimates based on 
the Sweet-Parker reconnection rate (Sweet 1958, 
Parker 1957). In terms of our simulations, we note 
that if reconnection is slow than we expect that 
magnetic fields will show a thin fibril structure in 
the saturated state, with a typical field reversal 
scale much smaller than the curvature scale. We 
see from the structure function that this is not 
the case. Apparently, in our simulations at least, 
reconnection is fast enough to relax tangled mag- 
netic field structures in no more than an eddy turn 
over time. Either the suppression of reconnection 
appears suddenly at higher resolutions, or theo- 
retical arguments suggesting fast reconnection in 
highly conducting turbulent fluids (Lazarian and 
Vishniac 1999) are correct. If the latter interpre- 
tation holds up, then arguments suggesting the 
suppression of the a effect (e.g. Vainshtein and 
Cattaneo 1992) at high Reynolds numbers are un- 
likely to prove correct. 

In what follows, we will give a rough explana- 
tion for the dependence of field quantities on v 
(=rj). Readers are advised that this analysis might 
fail in the limit of vanishing diffusivity. The reason 
is that run 144A shows slight deviations from the 
extrapolation results (see V 2 in Fig. 4a and T b \ V2 
in Fig. 6a). However, in most cases, the deviation 
is small. 

The kinetic energy depends on both viscosity 
and magnetic energy. The result shown in Fig. 
4a that the kinetic energy is almost independent 
of v = rj is actually the result of a rough cancella- 
tion between these two competing effects. Suppose 
that magnetic field is turned off and, therefore, 
the turbulence is ordinary hydrodynamic turbu- 
lence. Let us also suppose that the viscosity is 
zero. Then, at the energy injection scale, we will 
have 

V 3 ~ e. (22) 

This is a well known relation in Kolmogorov phe- 
nomenology wherein the energy injection rate e is 
equal to the energy cascade rate at the largest en- 
ergy containing scale, which is proportional to the 
energy contained in the largest eddies (~ V 2 ) di- 
vided by one eddy turnover time (~ L/V). When 
viscosity is not extremely small, we need to add a 
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viscosity effect: 



V 3 + vV 2 



(23) 



where we omit constants of order unity in each 
term. The second term on the left-hand side 
is proportional to dissipation at the large scale. 
Since the second term on the left-hand side is 
much smaller than the first, we can assume that 
vV 2 oc v + 0{v 2 ). If we solve this relation for V 2 , 
we obtain 



V 2 ~ C f e 2 ' 3 - O, 



(24) 



where the C's are constants. Now if we turn on 
the magnetic fields and set the viscosity to be zero, 
then we have 



V 3 + B 2 V ~e. 



(25) 



Here we also omit constants of order unity in 
each term. The first term on the left-hand side 
is still proportional to V 3 , which is supported by 
the relation T v i V2 oc V 3 (Fig. 6a). The second 
term on the left-hand side comes from the rela- 
tion T b i V2 oc VB 2 (Fig. 6b). Since the second 
term is small compared with the first, we can write 
B 2 V oc B 2 + 0{B 4 ). Solving for V 2 , we have 



V 2 



C e e 2 ' 3 - C B B 2 



(26) 



If we combine the two relations, we have 

V 2 ~ C € e 2/3 - C v v - C B B 2 , (27) 

where we ignore the interactions between the v- 
effect and the _B 2 -cffect. If any, such interactions 
will be weak because both the ^-effect and the B 2 - 
effect are weak. This equation implies that the 
quantity V 2 + CbB 2 depends only on v. 

To complete the above analysis, we need an ex- 
pression for B 2 . To this end, let us adopt an en- 
ergy cascade model analogous to the ordinary hy- 
drodynamical case: The energy transfer rate is 
equal to the magnetic energy cascade rate at the 
large scale. The energy transfer rate is given in 
equation (15). The magnetic energy cascade rate 
is B 2 /(L/V). Therefore, we have 



(V - B)B 2 ~ B 2 V + 7]B 2 



(28) 



where the second term on the right-hand side is 
the large scale dissipation term. Again, we omit 



all order unity constants. If we solve the equation 
for B, we obtain 



B = C V V - C n n, 



(29) 



where the C's are constants. In fact, our calcula- 
tion shows that 



5 -0.6V 



-26.7±2.1, 



(30) 



for Runs 128A, 96A, 72A, 64A, 64B, 64C, and 
64D. This implies that C v = 0.6 and C v ~ 26. 
This relation slightly underestimates the magnetic 
field strength for Run 144A. (It gives B — 0.6V, 
as opposed to the observed value of 0.8V.) If we 
examine Fig. 7 we see that this discrepancy arises 
from the difference between the rms magnetic field 
strength, and the magnetic field strength on the 
largest eddy scale. Comparing the hyperviscosity 
run (144A) with the other cases, we see that the 
former has a substantially greater fraction of its 
magnetic energy in small scale structures. If we 
interpret B in the preceding equation as the aver- 
age obtained by smoothing on the large eddy scale, 
which is consistent with our derivation, then the 
hyperviscosity run is in agreement with equation 
(30) . In this sense, the values of B and B 2 present 
in this discussion section should be regarded as 
lower limits for the case of vanishing diffusivity. 
Alternatively, one might regard the following dis- 
cussion as valid for cases of moderate magnetic 
Reynolds numbers and that this analysis sets lim- 
its for the case of zero diffusivity. For B 2 , we have 



B 2 ~ ClV 2 - 2C v C r yn. 



(31) 



Here we omitted C 2 n 2 on the right-hand side. 

If we combine expressions for V 2 and B 2 , we 
have 



B 



V 2 ~{C e e 2 l 3 -C v v + C v r,)/{l- 
(C 2 C t e 2 ' 3 - C 2 v C v v - C v C n Vn)/{\ 



■G B Cl), (32) 
-C B C 2 ). (33) 



Note that, even though we treat v and n sepa- 
rately, there is no guarantee that the above re- 
lations are true for non-unity Prandtl numbers. 
Since v = n and the second and third terms on 
the right-hand side of equation (32) have different 
signs, V 2 depends on v (=77) weakly. We expect 
Cb and Cy have similar values because they are 
describing similar things - dissipation. The ex- 
pression for B 2 is not absolutely correct because 
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the omitted term in equation (31) becomes non- 
negligible when r\ becomes large. 

When v = r\ = 0, we have 

V 2 ~ae 2/3 /(l + Ci?Ct>), (34) 
B 2 ~C 2 C £ e 2 / 3 /(l + C B C 2 ). (35) 

These give a lower limit for B 2 and an upper limit 
for V 2 . 

This work was partially supported by Na- 
tional Computational Science Alliance under 
CTS980010N and utilized the NCSA SGI/CRAY 
Origin2000. 
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Table 1 
Results of Simulations 



Run 


N 3 


V = 7] 








B 2 


e 


D M 


T i a 

V\V2 


T b\v 2 h 




(fl,f2) d 


REF1 


96 3 


.0043 


io- 


-4 


.773 


.169 


.167 


.088 


.0815 


.0353 




(500,560) 


REF2 


64 3 


.015 






0.85 




.163 










(300,750) 


REF3 


96 3 


.0043 






1.053 




.160 










(50,165) 


144A 


144 3 


hyper 


io- 


■7 


.649 


.420 


.161 




.0519 


.0758 


8.4 


(60,240) 


128A 


128 3 


.003 


io- 


■7 


.755 


.200 


.166 


.099 


.0794 


.0421 


7.1 


(100,210) 


96A 


96 3 


.0043 


io- 


■7 


.761 


.170 


.166 


.088 


.0807 


.0364 


6.7 


(200,500) 


72A 


72 3 


.0064 


io- 


■7 


.792 


.122 


.165 


.070 


.0870 


.0276 


6.8 


(300,800) 


64A 


64 3 


.0074 


io- 


■7 


.786 


.113 


.166 


.064 


.0835 


.0257 


6.6 


(300,800) 


64B 


64 3 


.01 


io- 


■7 


.794 


.0804 


.165 


.048 


.0860 


.0019 


6.6 


(300,750) 


64C 


64 3 


.015 


io- 


■7 


.808 


.0201 


.164 


.013 


.0903 


.0046 


8.0 


(300,800) 


64D 


64 3 


.02 


io- 


■7 


.774 


< io- 4 


.164 


- 10" 5 








(300,750) 



T V \V2 — -^t)|t>fc(^ 2) 

b T %2 =T %fc (fc = 2) 

c Equipartition wavelength 

d Time interval used for averaging physical quantities 
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Fig. 1. — Test of the code. Time evolution of Tu v , T v \ b , and magnetic dissipation. Tu v and — T v i b coincide 
exactly. Magnetic dissipation shows a time delay, which is the turbulent diffusion time scale. Run REF1. 

Fig. 2.— Time evolution of kinetic and magnetic energy The results at three different resolutions are shown: 
dotted, 144 3 (Run 144A); solid, 96 3 (Run 96 A); dashed, 64 3 (Run 64A). For 64A and 96A, the kinetic energy 
has very similar values. The magnetic energy depends on v (=77). In the case of hyperviscosity (Run 144A), 
the magnetic energy is more than half of the kinetic energy. 

Fig. 3. — Time evolution of kT b \ Vk {k) and E M (k). Run 144A 

Fig. 4. Dependence of physical quantities on v (= rj). Top: The kinetic energy density is nearly inde- 
pendent of the physical viscosity. The magnetic energy density grows as v (— 77) decreases. Bottom: This 
relation implies D M ~ (V - B/1.6)B 2 . 

Fig. 5. — Energy transfer and flux spectra for Run 144A. Top: Energy transfer spectra. The spectrum 
of T b \ Vk (dotted line) shows that magnetic fields gain energy from large scale eddies, bottom: Energy flux 
spectra. The spectrum of n b | tJfc (dotted line) rises rapidly at small k and reaches the value of e at k ~ 10. 

Fig. 6. — Evidences that large-scale turbulence is similar to ordinary hydrodynamic turbulence. Top: 
T v \ Vk (k = 2)/V 3 is nearly constant for nonzero physical viscosity, which implies that T v \ Vk (k = 2) <~ V 3 . 
In the case of hyperviscosity, this relation may not hold true. Bottom: T b \ Vh {k = 2)/(B 2 V). This relation 
implies stretching of magnetic field lines without significant back reaction is responsible for the nonlinear 
energy transfer from velocity components at k ~ 2 to magnetic fields. 

Fig. 7. — Energy spectra. The results at three different resolutions are shown: dotted, 144 3 (Run 144A); 
solid, 96 3 (Run 96 A); dashed, 64 3 (Run 64A). Both axes are drawn in logarithmic scale. 

Fig. 8. — Second order structure functions for the magnetic and velocity fields. The horizontal axis is in the 
direction of the local magnetic field, defined by pairwise averaging. The velocity fields show correlations at 
all scales, with a smooth transition from isotropy on large scales to anisotropy on small scales. The magnetic 
field shows similar behavior on small scales, below the equipartition scale, but has almost no correlation on 
large scales. 
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